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A precise measurement of the cross section of the process e + e~ — > 7r + 7r~(7) from threshold to 
an energy of 3 GeV is obtained with the initial state radiation (ISR) method using 232 fbr 1 of 
data collected with the BABAR detector at e + e~ center-of-mass energies near 10.6 GeV. The ISR 
luminosity is determined from a study of the leptonic process e + e~ — > /i + pi~7(7). The leading- 
order hadronic contribution to the muon magnetic anomaly calculated using the mv cross section 
measured from threshold to 1.8 GeV is (514.1 ± 2.2(stat) ± 3.1(syst)) x 10~ 10 . 

PACS numbers: 13.40Em, 13.60.Hb, 13.66. Be, 13.66.Jn 



Measurements of the e + e — > hadrons cross section are necessary to evaluate dispersion integrals for calculations 
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of hadronic vacuum polarization (VP). Of particular in- 
terest is the contribution a!^ ad to the muon magnetic mo- 
ment anomaly a^, which requires data in a region dom- 
inated by the process e + e~ — > 7r + 7r~(7). Comparison 
of the theoretical and measured [l| values of a M shows a 
discrepancy of about 3a when current e + e~ data 0,0, HI 
are used, possibly hinting at new physics. An approach 
using t decay data corrected for isospin-breaking, leads 
to a smaller difference 

The results on tttt production reported in this Letter 
are obtained with the ISR method [6( using e + e _ annihi- 
lation events collected at a center-of-mass (CM) energy 
y/s near 10.58 GeV. The cross section for e + e~ — > X 
at the reduced energy yfs' = mx, where X can be any 
final state, is deduced from a measurement of the radia- 
tive process e + e~ — > Xj where the photon is emitted by 
the e+ or e~; s' = s(l - 2E*/y/s), where £* is the CM 
energy of the ISR photon. In this analysis, v / s' ranges 
from threshold to 3 GeV. Two-body ISR processes with 
X = n + Tr~("/) and X — fi + fj," (7) are measured, where 
the ISR photon is detected at large angle and the charged 
particle pair can be accompanied by a final state radia- 
tion (FSR) photon. Obtaining the nn cross section from 
the ratio of pion to muon yield reduces significantly the 
systematic uncertainty. The measured muon cross sec- 
tion is compared to the QED prediction, and this cross 
check of the analysis is termed the QED test. 

The \fs' spectrum of e + e~ — > Xj events is related to 
the cross section for the process e + e _ — > X through 
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where 6x7 is the detection efficiency (acceptance) deter- 
mined by simulation with corrections obtained from data, 
and a x is the bare cross section (excluding VP). The 
measurement of a®^^ uses the effective ISR luminosity 

dLjg R / dy/s' provided by the measured mass spectrum of 
^7(7) events following Eq.(fT]) in which cr x (y/s') is the 
l«/i(7) bare cross section computed with QED [7]. For 
the QED test, the measurement of er° , ■, uses the effec- 

v * ' MH7) 

tive ISR luminosity definition as a product of the e + e _ 
integrated luminosity (L ee ), the radiator function 0], the 
ratio of detection efficiencies for the ISR photon in data 
and simulation (not included in 6x7), and the VP correc- 
tion (a(s') I a(0)) 2 . The radiator function, determined by 
the simulation, is the probability to radiate one or several 
ISR photons so that the produced final state X (exclud- 
ing ISR photons) has mass \fs' . 

This analysis is based on 232 flV 1 of data recorded 
with the BaBAR detector [§] at the PEP-II asymmetric- 
energy e + e~ storage rings. Charged-particle tracks are 
measured with a five-layer double-sided silicon vertex 
tracker (SVT) together with a 40-layer drift chamber 
(DCH) inside a 1.5 T superconducting solenoid magnet. 
The energy and direction of photons are measured in the 



CsI(Tl) electromagnetic calorimeter (EMC). Charged- 
particle identification (PID) uses ionization loss dE/dx 
in the SVT and DCH, the Cherenkov radiation detected 
in a ring- imaging device (DIRC), and the shower deposit 
in the EMC (E ca i) and in the instrumented flux return 
(IFR) of the magnet. 

Signal and background ISR processes are simulated 
with Monte Carlo (MC) event generators based on 
Ref. 0| . Additional ISR photons are generated with 
the structure function method [l(|, and additional FSR 
photons with PHOTOS 11]. Background events from 
e + e~ — > qq (q = u, d, s, c) are generated with JET- 
SET [13]. The response of the BABAR detector is sim- 
ulated with GEANT4 

Two-body ISR events are selected by requiring a pho- 
ton with E* > 3 GeV and laboratory polar angle in the 
range 0.35 — 2.4 rad, and exactly two tracks of opposite 
charge, each with momentum p > 1 GeV/c and within 
the angular range 0.40 — 2.45 rad. If several photons are 
detected, the ISR photon is chosen to be that with the 
highest E*. The charged-particle tracks, required to have 
at least 15 hits in the DCH, must originate within 5 mm 
of the collision axis and extrapolate to DIRC and IFR 
active areas which exclude low-efficiency regions. An ad- 
ditional criterion based on a combination of E ca i and 
dE/dx reduces electron contamination. 

Acceptance and mass-dependent efficiencies for trig- 
ger, reconstruction, PID, and event selection are com- 
puted using the simulation. The ratios of data and MC 
efficiencies have been determined from specific studies, as 
described below, and are applied as mass-dependent cor- 
rections to the MC efficiency. They amount to at most a 
few percent and are known to a few permil level or better. 

Tracking and PID efficiencies are determined taking 
advantage of pair production. For tracking studies, two- 
prong ISR candidates are selected on the basis of the 
ISR photon and one track. A kinematic fit yields the 
expected parameters of the second track. The unbiased 
sample of candidate second tracks is used to measure 
track reconstruction efficiency. The maximum correlated 
two-track loss induced by track overlap in the DCH is 
0.6% for pions and 0.3% for muons. 

Tracks are assigned uniquely to a complete set of PID 
classes using a combination of cut-based and likelihood 
selectors. The '/it' class is addressed first by making use of 
track IFR penetration and hit spread distribution, and 
of the E ca i value. Tracks failing the '/i' identification 
are labeled as 'e' if they satisfy E ca i/p > 0.8. The 'X' 
class is determined using DIRC information and dE/dx. 
Remaining tracks are labeled as V. A tighter selection 
called '71V is applied in mass regions where background 
dominates or to create a pure pion test sample. 

Efficiencies for PID are measured from pure samples 
of muon, pion, and kaon pairs obtained from xxj events 
where one track is selected as 'ir^, or 'K' and the 
other is used to probe the PID algorithm. The efficien- 
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cies are stored, according to momentum and position in 
the IFR or the DIRC. Typical efficiency for y is 90%, 
with 10% mis-ID as l ir'. The V efficiency is strongly 
momentum-dependent because of mis-ID as i K' (1% at 
lGeV/c, reaching 20% at 6GeV/c), as y (5-6%), or as 
'e' (2%). Correlations between the PID efficiencies, due 
to track overlap, have been observed and parametrized. 
They are largest for muons where the correlated PID loss 
reaches 1.3% of the events below 1 GeV/c 2 . It is impor- 
tant to control this effect, since it affects the irw and pp 
samples in an anti-correlated way. 

To obtain the spectra Njj of produced particle pairs of 
true type j, a set of three linear relations must be solved. 
They involve the Njj, the measured mass distributions 
for each Hi '-identified final state, and the probabilities 
£ °4i> {hi — Mi n or K) which represent the product of the 
measured efficiencies for each track of true type j to be 
identified as H\ corrected by correlation factors. 

A contribution (< 10~ 3 ) to N V7T from ppj is esti- 
mated from MC and subtracted after reweighting the 
rate to agree with the BABAR measurement [14j. Multi- 
hadronic background from e + e~ — * qq comes from low- 
multiplicity events in which an energetic 7 originating 
from a ir° is mistaken as the ISR photon candidate. To 
normalize this rate from JETSET, the 7r° yield obtained 
by pairing the ISR photon with other photons in the 
event is compared in data and MC; JETSET overesti- 
mates this background by a factor 1.3. Multi-hadronic 
ISR backgrounds are dominated by e + e~ — ► 7r + 7r~7r°7 
and e + e~ — * 7r + 7r~27r°7 contributions. An approach 
similar to that for qq is followed to calibrate the back- 
ground level from the 37r ISR process, using ui and cj> 
signals. The ratio of data to MC yield is found to be 
0.99 ± 0.04. The MC estimate for the 27r27r°7 process 
is used and assigned a 10% systematic uncertainty. A 
residual radiative Bhabha background is identifiable only 
near threshold and at large mass, where the pion signal 
vanishes. Its magnitude is estimated from the helicity 
angle distribution in the tttt CM frame at low energy and 
its energy dependence obtained from a control sample of 
radiative Bhabha events. It is assigned a 100% system- 
atic uncertainty. To suppress the contribution from the 
e + e~ — > 77 process with a photon conversion, which af- 
fects the spectrum at threshold, the vertex of the two 
tracks is required to be closer than 5 mm to the collision 
axis in the transverse plane. This criterion is applied only 
to events in the p tails, defined to lie outside the central 
region 0.5 < m vv < 1.0GeV/c 2 . Background contribu- 
tions to the N^n spectrum are negligible. 

Each event is subjected to two kinematic fits to the 
e + e _ — > Xj hypothesis, where X includes one addi- 
tional photon, detected or not. Both fits use the ISR 
photon direction and the parameters and covariance ma- 
trix of each charged-particle track. The energy of the 
ISR photon is not used, as it has little impact for the 
relatively low CM energies involved. The two-constraint 



(2C) 'ISR' fit allows an undetected photon collinear with 
the collision axis, while the 3C 'FSR' fit uses any photon 
with E 1 > 25 MeV. When more than one such photon 
is present, the best 'FSR' fit is retained. An event with 
no extra photon is characterized only by its Xisr va l ue - 
Most events have small x 2 values for both fits; an event 
with only a small Xisr {X%sr) indicates the presence 
of additional ISR (FSR) radiation. Events where both 
fits have large % 2 values result from track or ISR pho- 
ton resolution effects, the presence of additional radiated 
photons, or multi-hadronic background. To accommo- 
date the expected background levels, different criteria in 
the (XisR'Xfsr) pl ane are applied depending on the m n „ 
mass regions. For the central p region, a loose 2D contour 
has been optimized to remove the main background area 
while maintaining control of the associated systematic 
uncertainties. The same procedure is used in the ppj 
analysis in spite of the very small background. In the p 
tails, a tighter x 2 selection is imposed to reduce the larger 
background. Samples of 529320 pion and 445631 muon 
events are selected in the mass range below 3GeV/c 2 , 
where the m nw (m^) mass is calculated from the best 
'ISR' or 'FSR' fit. 

The computed acceptance and the \ 2 selection effi- 
ciency depend on the description of radiative effects in 
the generator. The FSR rate is measured from events 
that satisfy the 'FSR' fit, with an additional photon 
(E 7 > 0.2 GeV) within 20° of either track. The ex- 
cess in data relative to the generator prediction using 
PHOTOS [ll| is (-4 ± 6)% of total FSR for muons, 
and (21 ± 5)% for pions. This difference results in a 
(6 ± 2) x 10 -4 correction. More significant differences 
are found between data and the generator for additional 
ISR photons, since the latter uses a collinear approxima- 
tion and an energy cut-off for very hard photons. In- 
duced kinematical effects have been studied using the 
next-to-leading order (NLO) PHOKHARA generator [l6| 
at four-vector level with fast simulation. Differences in 
acceptance occur at the few percent level, and these yield 
corrections to the QED test. In contrast, since radia- 
tion from the initial state is common to the pion and 
muon channels, the mr^) cross section, obtained from 
the irir/pp ratio, is affected and corrected only at a few 
permil level. Additional ISR effects on the \ 2 selection 
efficiencies factorize in both processes and cancel in the 
ratio. The \ 2 selection efficiency determined from muon 
data applies to pions, after correction for the effect of 
secondary interactions and the n/fi difference for addi- 
tional FSR. Therefore the measurement of the pion cross 
section is to a large extent insensitive to the description 
of NLO effects in the generator. 

The QED test involves two additional factors, both of 
which cancel in the tttt/ pp ratio: L ee and the ISR photon 
efficiency, which is measured using a ppj sample selected 
only on the basis of the two muon tracks. The QED test 
is expressed as the ratio of data to the simulated spec- 
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trum, after the latter is corrected using data for all known 
detector and reconstruction differences. The generator is 
also corrected for its known NLO deficiencies using the 
comparison to PHOKHARA. The ratio is consistent with 
unity from threshold to 3GeV/c 2 , (Fig. Q] (a)). A fit to a 
constant value yields (x 2 / n df = 55.4/54; n c if=number of 
degrees of freedom) 



S 105 
o 



<7 



data 
W7(7) 



NLO QED 



1 = (40 ± 20 ± 55 ± 94) x 10" 



(2) 



where the errors are statistical, systematic from this anal- 
ysis, and systematic from L ee , respectively. The QED 
test is thus satisfied within an overall accuracy of 1.1%. 

To correct for resolution and FSR effects, an unfold- 
ing of the background-subtracted and efficiency-corrected 
TOtttt distribution is performed. A separate mass-transfer 
matrix is created using simulation for the p central and 
tail regions; this provides the probability that an event 
generated in a y/s' interval i is reconstructed in a m„ 
interval j. The matrix is corrected using data to account 
for the larger rate of events with poorer mass resolution. 
Performance and robustness of the unfolding method [l7| 
have been assessed using test models. For the 2-MeV 
intervals, the significant elements of the resulting covari- 
ance matrix lie near the diagonal over a typical range of 
6 — 8 MeV, which corresponds to the energy resolution. 

The results for the e + e~ 
tion including FSR, a% 

(b). Prominent features are the dominant p resonance, 
the abrupt drop at 0.78 GeV due to p — u interference, a 
clear dip at 1.6 GeV resulting from higher p state inter- 
ference, and additional structure near 2.2 GeV. System- 
atic uncertainties are estimated from the precision of the 
data-MC comparisons and from the measurement proce- 
dures used for the various efficiencies. They are reported 



TABLE I: Relative systematic uncertainties (in 10 ) on the 
e + e~ — > 7r + 7r~(7) cross section by y 7 ? intervals (in GeV) up 
to 1.2 GeV. The statistical part of the efficiency uncertainties 
is included in the total statistical uncertainty in each interval. 



7T + 7r (7) bare cross sec- 
(7) (v / s 7 ), are given in Fig. CD 



Source of 


CM Energy Interval (GeV) 


Uncertainty 


0.3-0.4 


0.4-0.5 


0.5-0.6 


0.6-0.9 


0.9-1.2 


trigger/ filter 


5.3 


2.7 


1.9 


1.0 


0.5 


tracking 


3.8 


2.1 


2.1 


1.1 


1.7 


7T-ID 


10.1 


2.5 


6.2 


2.4 


4.2 


background 


3.5 


4.3 


5.2 


1.0 


3.0 


acceptance 


1.6 


1.6 


1.0 


1.0 


1.6 


kinematic fit {\ 2 ) 


0.9 


0.9 


0.3 


0.3 


0.9 


correlated pp ID loss 


3.0 


2.0 


3.0 


1.3 


2.0 


rnv/pp non-cancel. 


2.7 


1.4 


1.6 


1.1 


1.3 


unfolding 


1.0 


2.7 


2.7 


1.0 


1.3 


ISR luminosity (pp) 


3.4 


3.4 


3.4 


3.4 


3.4 


total uncertainty 


13.8 


8.1 


10.2 


5.0 


6.5 



I - { f i i * 

nqs U (a) 




2 2.5 3 

nJs 7 [GeV] 

FIG. 1: (a) The ratio of the measured cross section for 
e + e~ -> p + p~l{l) to the NLO QED prediction. The 
band represents Eq. ([2]). (b) The measured cross section for 
e + e~ —> tx + tx~ (7) from 0.3 to 3 GeV. (c) Enlarged view of the 
p region in energy intervals of 2 MeV. The errors are from the 
combined diagonal elements of the statistical and systematic 
covariance matrices. 



in Table U for 0.3 < W < 1.2 GeV. Although larger 
outside this range, the systematic uncertainties do not 
exceed statistical errors over the full spectrum for the 
chosen energy intervals. 

The lowest-order contribution of the 7777(7) intermedi- 
ate state to the muon magnetic anomaly is given by 



r(~f),LO __ 



OO 

4m?. 



(*') , (3) 



where K(s') is a known kernel [19[. The integration uses 
the measured cross section and the errors are computed 
using the full statistical and systematic covariance ma- 
trices. The systematic uncertainties for each source are 
taken to be fully correlated over all mass regions. The 
integrated result from threshold to 1.8 GeV is 



o™ (7) ' i0 = (514.1 ± 2.2 ± 3.1) x 10" 



10 



(4) 



where the errors are statistical and systematic. This 
value is larger than that from a combination of previ- 
ous e + e" data [B| (503.5 ±3.5), but is in good agreement 
with the updated value from t decay [B| (515.2 ± 3.4). 

In summary, the cross section for the process e + e~ — > 
7r + 7r~ (7) has been measured in the energy range from 0.3 
to 3 GeV, using the ISR method. The result for the 7T7t 
hadronic contribution to has a precision comparable 
to that of the combined value from existing e + e~ experi- 
ments. However, the BABAR central value is larger, which 
reduces the deviation of the direct a M measurement from 
the Standard Model prediction. 
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